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Abstract 



Extending our previous work on SU(3), we construct spherically symmetric 
BPS saturated regular configurations of M = 4 SU(N) supersymmetric Yang-Mills 
theory preserving 1/4 supersymmetry, and investigate their features. We also give 
exact solutions in the case some of the free parameters of the general solutions take 
certain values. These field theory BPS states correspond to the string theory BPS 
states of multi-pronged strings connecting N different D3-branes by regarding the 
M = 4 supersymmetric Yang-Mills theory as an effective field theory on parallel 
D3-branes. We compare our solutions with multi-pronged strings in string picture. 
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1 Introduction 



In recent developments of non-perturbative string theory, the existence of many kinds of 
extended BPS objects has been realized. These BPS states are related to each other by several 
duality transformations, which are believed to be the true symmetries of the non-perturbative 
string theory. By studying a supersymmetric gauge theory as an effective field theory on an 
appropriate brane configuration, the existence of some string theory BPS states on the brane 
configuration implies the existence of the corresponding BPS states in the field theory. Thus 
showing explicitly the existence of such corresponding BPS states by using traditional field 
theory techniques would be a good test to the non-perturbative implications of string theory 
and would deepen the understanding of the field theory itself. 

A typical example of such interplay is given by the 3+1 dimensional Af = 4 SU(N) super- 
symmetric Yang-Mills (SYM) theory broken spontaneously to ^(l)^ -1 . This theory can be 
studied as an effective field theory on N parallel D3-branes [p], 3- The invariance of D3-branes 
under the SL(2, Z) duality transformation of type IIB string theory implies the SL(2, Z) du- 
ality symmetry of the M = 4 SYM theory 0, |[]. The (±1,0) strings stretching between 
different D3-branes preserve 1/2 supersymmetry of the D3-brane world volume and appear 
as the massive W-bosons of the broken gauge symmetries in the field theory. The SL(2, Z) 
duality symmetry implies the existence of other BPS states of the field theory corresponding 
to the general (p, q) strings with relatively prime integers p and q. The states with (p, ±1) 
are monopoles and dyons, the field configurations of which are explicitly known in the form of 
the Prasad- Sommerfield solution ||. The existence of the states with q = 2 was shown by the 
quantization of the collective modes of the two monopole solutions ||], while the existence of 
the states with q > 2 was discussed in 0. 

The BPS objects we discuss in this paper are multi-pronged strings in type IIB string 
theory and their counterparts in M = 4 SYM theory. It was conjectured originally in [|7| 
that a three-pronged string would exist as a stable state under both the conditions that the 
forces from the three strings should balance and that the two-from charges of the three strings 
should conserve at the junction. The BPS feature of such a junction has been shown by 



several authors in both the string picture M, |J and the M-theory picture fllO |. The relations 



to the field theory BPS states have been investigated in several contexts [11, 12, 13, 14]. In 



12|, a three-pronged string connecting different parallel D3-branes was investigated. This 



configuration preserves the 1/4 of the D3-brane world volume supersymmetry, and hence it 
is expected that there exist corresponding 1/4 BPS states in the M = 4 SU(N) SYM theory 
with N > 2. The author of [12| obtained an evidence for the existence, showing that the 
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mass of the three-pronged string agrees with the mass of the corresponding field theory state 
under the assumption of its BPS saturation. But there the configuration itself in the field 
theory was not obtained. Such a 1/4 BPS state was shown to have non-parallel electric and 



magnetic charges (JT2J, jig] , and this fact makes it a non-trivial issue to solve the equations of 
BPS saturation. 

In our previous work fl6f , we constructed explicitly the corresponding BPS states in the 
M = 4 SU(3) SYM theory. We solved explicitly the Bogomol'nyi equations in the case of 
non-parallel electric and magnetic charges. We also obtained exact solutions in the case some 
of the free parameters are fixed to certain values. The asymptotic behavior of our solutions 
agrees with the string picture, that is, three straight strings meet at one point and the forces 
balance at the junction. On the other hand, we pointed out that there are some differences 
in our solutions from the string picture, and the behavior of our solutions near the junction 
does not seem to respect the SL(2, Z) duality symmetry. This fact may be interpreted as our 
technical preference that we treat the SYM theory classically, i.e. electrically. 

In this paper, we continue our investigations of string junctions from field theoretical 
viewpoints, and extend our analysis to multi-pronged strings. The main part of this paper 
consists of the explicit construction of the 1/4 BPS states in M = 4 SU(N) SYM theory 
(N > 2). We solve the Bogomol'nyi equations under the assumption of spherical symmetry 
in the D3-brane world volume. The BPS solutions we obtain correspond to the string theory 
BPS states of an iV-pronged string connecting N different D3-branes. Such multi-pronged 
strings have been recently studied in the string picture in |TJ . 



This paper is organized as follows. In section |2|, we explain our strategy to solve the Bogo- 
mol'nyi equations. We investigate some properties of our solutions, and give exact solutions 
in the case some of the free parameters of the general solutions take certain values. In section 
H we show in detail the features of our exact solutions for SU(4) and SU(5) as examples. 
Various interesting properties of the solutions are also discussed there. In section f|, we calcu- 
late the long-range force between two dyons in the SYM theory, and find the consistency with 
the idea that a three-pronged string can be produced as a bound state of two strings. Section 
D contains summary and discussions. In appendix [E], we show that our solutions satisfy the 
equations of motion of the non-Abelian Born-Infeld action, and the other appendices contain 
several formulas used in the text. 
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2 BPS saturated SYM configurations 



2.1 BPS bound and Bogomol'nyi equations 



As mentioned in sec. [I], we would like to generalize the arguments of ref. [16] for the 3 + 1 



dimensional M = 4 SYM with gauge group SU (3) to the SU (N) case and construct the BPS 
saturated solutions which correspond to the string theory BPS states of multi-pronged strings 
connecting N D3-branes. Omitting the fermionic terms, the energy of our SU(N) SYM system 
reads, 

U = J d 3 x^ Tr {(B,) 2 + {Si) 2 + (AY) 2 + (AY) 2 + (D X) 2 + (D Y) 2 - [X, Y] 2 } , (2.1) 

where Si = and Bi = Xj^e^F^ are the electric and magnetic fields, and the covariant 
derivative is defined by D^X = d^X — i [A^, X]. We have put the Yang-Mills coupling constant 
equal to one, and shall consider the case of vanishing vacuum theta angle. In (|2.1|) , we have 
kept only two adjoint scalars, X and Y, among the original six which should describe the 
transverse coordinates of the D3-branes. This is because we are interested in the multi- 
pronged BPS string states preserving 1/4 of the supersymmetry, and such string networks 
must necessarily lie on a two-dimensional plane. 

In this paper we assume the D-flatness [X, Y] = 0.* To derive the BPS saturation condition, 



we rewire the energy (|2.1|) as [15 



U = J d 3 x^ Tr {(Si cos 9-Bi sm9- DiX) 2 + (B l cos 9 + Si sin 9 - DiY) 2 

+ (A)X) 2 + (A)Y) 2 } + {Qx + My) cosfl + [Q Y - M x ) sm#, (2.2) 
where 6 is an arbitrary angle, and Qx,y and M X y are defined by 

Qx = J d 3 x Tr (SiDiX) = J dS { Tr , M x = J d 3 x Tr (^ AY) = J dS { Tr (BiX) . 

r— >oo r—*oo 

(2.3) 

Since each term in the integrand of ( p.2[ ) is non-negative and 9 is arbitrary, we obtain the BPS 
bound: 

U > Abps = ^(Qx + M Y ) 2 + (Q Y -M x ) 2 . (2.4) 
The lower bound in ( |2.4| ) is saturated when the following conditions hold: 

DiX = Si cos 9 - Bi sin 9, (2.5) 

* Due to the D-flatness, we can diagonalize X and Y simultaneously and hence consider definite D3-brane 
surfaces described by the eigenvalues of X and Y as we shall do in later subsections. The Bogomol'nyi equations 
without imposing the D-flatness are given in refs. p% pq] 
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with the angle 9 given by 



DiY = BiCosO + Si sine, (2.6) 
D X = D Y = 0, (2.7) 
[X,Y\ = 0, (2.8) 



. Qy-M x Qx + M y 

sm.0 = — , cosO = — . (2.9) 

Ebps -Ebps 



In addition to the four equations (|2.5| )- (|2.8| ), we have to impose the Gauss law, 

DiSi = 0, (2.10) 



since we used it in converting the volume integration of Qx,y ( P-3| ) into the surface one. Note 
that eq. ( |2.9| ) is an automatic consequence of the two equations (|2.5| ) and ( |2.6| ) and need not 
be imposed independently. In appendix [D] we show that the SYM configurations satisfying 
eqs. (|2.5|) - (|2.8|) preserve 1/4 supersymmetry. 



2.2 SYM solutions and the IIB picture of string networks 

In this subsection, we shall discuss the asymptotic (r — > oo) behavior of the solutions and 
its relation to the IIB picture of the string networks. Suppose that we have a static solution 
(A^r), X(r), Y(r)) to the equations (|2.5|) - (|2.8|) and (|2.10|) , and that their asymptotic forms 
are, after a suitable gauge transformation, given (locally) as follows: 

X ~ diag(x a ) + — diag(w a ), (2.11) 
2r 

Y ~diag(y a ) + ^-diag(t; a ), (2.12) 
2r 

T ' I 

Si |-diag(e a ), (2.13) 

r z 2 

B t ^^ldmg(g a ), (2.14) 
r z 2 

with r = r/r. For these asymptotic behaviors, we have 

Af N 

Qx = 27r e a Xa, M x = 2ti ^ g a x a , (2.15) 

a=l a=l 

and similarly for Qy and My. 

Now, this solution is interpreted as representing a configuration of iV D3-branes a = 
1, 2, . . . , iV at transverse coordinates (x a , y a ) from which a string with two-form (NS-NS and 
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D3-brane 

Figure 1: "Tube-like" configuration of D3-brane surface representing a string 



R-R) charges (e a ,g a ) are emerging in the direction (u a ,v a ). This is because the eigenvalues 
of the scalars (X, Y) are interpreted as the transverse coordinates of the D3-branes and the 
"tube-like" part of the D3-brane surface (corresponding to smaller r) can be regarded as a 
string |T7| (see fig. |T]). The string directions (u a , v a ) are not arbitrary but are related to (e a , g a ) 

byi 

(2-16) 



cosy 
sin 9 



— sin 6 

cos 6 



9a 



where 9 is the angle given by ( |2.9|) . Leaving the derivation of eq. fl2.16|) for our SYM solutions 
in sec. |2l|, we shall discuss the IIB picture of string networks deduced from fl2.16|) . 



In fact, we can show that, given (e a , g a ) and (x a , y a ) (a — 1, 2, . . . , N) satisfying eqs. (|2.9| ) 
and (|2.16| ), it is possible to draw "generalized" tree string networks made of 3-string junctions 
satisfying the following two properties: 



(A) The positions of the 3-string junctions in the network is consistently given in such a 
way that the direction of each component string A is parallel or anti-parallelQ to the 
vector (ua,va) which is related to the charge vector (e^,^) by eq. ( |2.16| ). Here, the 
component strings include the N strings emerging directly from the D3-branes as well as 
the internal strings connecting 3-string junctions, and the charges (eA,gA) are uniquely 
determined by the conservation condition at the junctions. 



t A string A connecting two 3-string junctions (or a D3-brane and a junction) a and (3 with coordinates 
{x a ,y a ) and {xp^p) is parallel (anti-parallel) to the vector {ua,va) defined to be directed from a to j3 if tA 
satisfying {x a ,y a ) + t A {u A , v A ) = {xp,Vp) is positive (negative). 
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(B) The sum of the masses of the component strings constituting the network coincides with 
the BPS bound E BPS ( p.4j) of the energy: 

£bps = E T ^4- ( 2 - 17 ) 

A 

Here, T4 = ±2ii^J e\ + g\ and £a are the tension and the length of the string A, re- 
spectively, and the sign factor of the tension Ta is positive (negative) if the string A is 
parallel (anti-parallel) to the vector (ua,Va)- 

In appendix [A|, we present a proof of the properties (A) and (B) valid for general N using 
a technique which reduces the problem to the simple N = 3 case. The string networks 
satisfying (A) and (B) are "generalized" ones in the sense that we are allowing strings with 
negative tension (see fig. |2|). As seen from the proof given in appendix [A], there are in general 




(a) (b) 

Figure 2: String networks corresponding to the same set of (e a ,g a ) and (x a ,y a ). The blobs 
represent D3-branes. The network (a) is physical, while (b) is unphysical since it contains a 
string with negative tension (dashed line). 

more than one possible string networks for a given set of (e a ,g a ) and (x a ,g a ). In order for it 
to be possible to identify a physical network consisting only of strings with positive tension, 
(x a ,y a ) cannot be arbitrary but have to satisfy conditions determined by the charges (e OJ g a ). 
Such a condition is given by eq. (|4.1j) for the N = 3 case. In sec. [3.1| , we analyze our exact 
SYM solutions satisfying the BPS saturation condition for the N = 4 case to find that for 
every parameter value of the solution there are corresponding physical networks. Though 
it is generally expected that, for (x a ,y a ) and (e a ,g a ) associated with a BPS saturated SYM 
configuration, we can always identify physical networks, we do not yet have a proof of it. 

We mentioned at the beginning of this subsection the interpretation of the eigenvalues 
of the scalars (X, Y) as the coordinates of the IIB strings emerging from the D3-branes. 
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However, as we saw in ref. [16 for the N = 3 case and shall find later in this paper for 
general N, such interpretation does not match the string network picture for finite r since the 
lines corresponding to the eigenvalues of the scalars are generally curved. Moreover, multiple 
3-string junctions in the networks with iV > 4 is impossible in this interpretation due to the 
analyticity of the scalars as functions of the four dimensional coordinates of the SYM. 



2.3 Spherical symmetry and monopole solutions 



Let us prepare a framework for constructing the BPS saturated SYM configurations. Since 
the angle 9 can be absorbed by a rotation in the (X, Y) plane, we shall consider the case of 
9 = 0. Then, the equations to be solved are 



DiX — Si, 



(2.18) 
(2.19) 



as well as ( |2.7|) , ( |2.8| ) and ( [2.10|) . Our strategy for the construction of the solutions is the same 
as in ref. |16| for the SU(3) case. First, we prepare an SU(N) monopole solution (A^r), Y(r)) 
to eq. Q2.19| ). Then, eq. fl2.18| ) is automatically satisfied by putting A (r) = —X(r), while eq. 
(|2~?1 ) holds due to eq. ( |2.8|) and the time-independence of our solution. Therefore, we have 
only to solve eq. (|2.10| ), i.e., 

DMX = 0, (2.20) 

under the D-flatness condition 



). As we did in ref. [[RJ], we adopt as the monopole solutions 
to eq. (|2.19|) those given in ref. |1| constructed using a general formalism for spherically 
symmetric solutions of ref. [|I9"| |. In the following, we shall present our problem in the formalism 
of ref. [[Lj| and explain the monopole solutions of ref. |18|] (see also ref. [^(]] ) . 

Our solutions (Ai,X, Y) are assumed to be spherically symmetric with respect to an an- 
gular momentum operator J, i.e., they satisfy 



[J u Aj] = ie ijk A k , [Ji, X] = [Jj, Y] = 0, 
The present J is the sum of the space and the gauge group rotations: 

J = L + T, 



(2.21) 



(2.22) 



where L = —ir x V is the generator of the space rotation and T is the maximal SU(2) 
embedding in SU (N) with T 3 = | diag (N — 1, N — 3, • • • , —N + 1). The monopole solution 



of |T8[ assumes the following form for the vector potential: 

A(r) = (M(r, f ) - T) x r/r, 



(2.23) 
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where the Lie algebra valued function Mj should satisfy the spherical symmetry condition, 



[J i ,M j 



ie ijk M k . 



(2.24) 



Various formulas are derived by using the expression V = rd/dr — {i/r)r x L for the space 
derivative as well as the spherical symmetry properties, eqs. ( |2.21|) and ( |2.24| ). We need in 
particular the following three: 



DY = rY' + - f x [M, Y] , 



r 

DiDiX = X" + -X'-\ (<% - hf 3 ) [M is [Mj, X]] 



(f x (f x M')) 



i i 



(2.25) 
(2.26) 
(2.27) 



where the prime denotes the differentiation d/dr. 

Due to the spherical symmetry it is sufficient to construct solutions on the positive z-axis. 
The monopole equation ( |2.19| ) and the Gauss law constraint fl2.20| ) on the z-axis are reduced 
to 
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r z Y' = -[M + ,M_]-T 3 , 



M' ± = T [M ± ,Y], 

X" + -X' - r^([M + , [M.,X\] + [M_, [M+,X\]) = 0, 



(2.28) 
(2.29) 

(2.30) 



with M± = Mi ± iMi- As seen from eqs. ( |2.18| ), ( |2.19| ) and (|2.25| ), the ^-component of the 
electric and magnetic fields on the positive z-axis are equal to the derivatives of the scalars: 



S z = X\ 



B z = Y'. 



(2.31) 



The matrix-valued functions M±(r) and Y{r) of ref. [TB| are given on the positive z-axis 



as 



(M + ) mr 

F = diag(F m (r)) 



$m,n— 1 ^mijj 
1 



Ml, 



2r 



diag(^ m (r) -* m _i(r)), # c 



N 



using a m and \l/ m , which are further expressed in terms of Q m {r) as 



,1/2 



/ 

a m( r ) = 7^— ( mm Qrn~lQm+l 
^m( r ) = ~ r (hlQm)' + ™7l = N — m), 



(2.32) 
(2.33) 

(2.34) 
(2.35) 
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with Qq — Q N — 1. The functions a m and M± satisfy eqs. Q2.28|) and (|2.29|) if Q m ( r ) is a 
solution to 

{Q'mf ~ QmQm = mm Qm+lQm-1, (2.36) 

for m — 1,2, . . . ,N — 1. In ref. [pi| they found the following Q m satisfying (|2.36| ): 

(2.37) 



Q m =i m j2 n e ~ 2yar n (2y6-2y a )- 1 



where 7 m is given by 



lr, 



(2.38) 



and the sum in ( |2.37| ) is over the distinct ways of dividing the integers {1, 2, ... , N} into 
two groups, D m with m elements and D m with m elements. 



In ( |2.37D , y a (a = 1, 2, . . . , N) are arbitrary parameters satisfying J2a=i Va — 0. They are 



nothing but y a appearing in the asymptotic expression (|2.12|) of the scalar Y if the condition 
Vi < Vi < • • • < Vn is satisfied. This is seen from the asymptotic behavior of hiQ m , 



In Q m = -2 ^2 Va r + const. + O (e~ 



(r — > oo 



(2.39) 



a=l 



and eq. ( |2.35| ). Using ( |2.31| ), we find that the asymptotic form of the magnetic field is B z 
—Tz/r 2 and hence the magnetic charges g a of the present solution are given by 



((/ a ) = (-JV + l,-JV + 3,...,JV-l). 



On the other hand, the behavior of Q m (r) of (|2.37f ) near the origin r = is 

Q m = r mm (l + q m r 2 + 0(r 3 )) (r ~ 0), 

where the coefficient q m of the sub-leading term is given byQ 

2mm 



(2.40) 



(2.41) 



N 



N(N 2 - 1) ^ 

The leading behavior r mm of Q m is consistent with the regularity of the solution at r = 0. 



(2.42) 



* One way to derive eq. ( 2.42[ ) is as follows. Eq. ( 2.36 ) implies that q m satisfies the recursion relation 
(2mm- 2)q rn = mm(q m+ i —g m -i), whose solution is given by q m oc mm (incidentally, this recursion equation 
is the p — 1 case of eq. ( 2.48| ), which is solved in appendix [§]). Taking into account the initial condition 
9i = jv(jv+i) E^Li vl obtained from ( p.37[ ), we get ( $1% ). 
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Since we have adopted the diagonal form (|2.33| ) for the scalar Y, the D-flatness condition 
(|2.8|) implies that the other scalar X is also diagonal. The diagonal form of X is also a 
consequence of the spherical symmetry [J3,X] = fl2.21|) implying [T 3 ,X] = on the z-axis. 
Therefore, we express X in terms of TV — 1 functions $ m (m = 1, 2, . . . , N — 1) as 

X = diag (X m ) = ±- diag ($ m - $ m _0 , $ = $ w = 0. (2.43) 
2r 

Then, eq. ( |2.3U| ) for X is reduced to 

< - ^ (2$ m - $ m+ i - $ m -i) = 0, m = 1, . . . , N - 1. (2.44) 

This differential equation can be solved numerically and for a certain special values of (y a ) 
analytically. In the next subsection, we shall present an exact solution which is a generalization 



of the exact solution for the N = 3 case given in ref. |16 |. 

Before closing this subsection, we derive the relation (|2.16|) for our SYM solutions. This 
is an immediate consequence of ( |2.31| ) for E z and B z and the asymptotic expressions ( [2.1 1| ) 
- ( gTTg ). They lead to eq. ( \LWj ) with 9 = 0, (u a ,v a ) = -(e a ,g a ). (Note that, since M±(r) 
decays exponentially as r — > oo as seen from eqs. ( |2.34| ) and ( |2.39| ), so do the other components 
S X) y and B x>y .) Eq. (|2.16| ) in the case of non- vanishing 9 is obtained by a rotation of (X, Y) 
and hence (u a , v a ). 



2.4 Exact solutions for X 



Similarly to the case of SU (3) fll6| , it is possible to construct an exact solution to the equations 
( |2.44j ) for the following special values of (y a ): 

(yi, y 2 ,...,y N ) = C (-N+ 1, -JV+3, ...,N-1), (2.45) 

where C is a positive constant. For this (y a ), Q m is given explicitly by0 

^ / sinh Cr\ mm 

Q m = f — — J , (2.46) 



and eq. ( ^.44| ) becomes 



TflTfl 

$ m - —TT- (2*m " " $m-l) = 0, m = 1, • • • , N- 1, (2.47) 

sinh r 



§ The easiest way to derive ( |2.46| ) is to first obtain Qi from eq. ( p. 37 ) and then use the recursion relation 
( 2.36| ) to get Q m with higher m. 
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where we have put C — 1 since C can be absorbed into the rescaling of r. 

To rewrite the differential equation (|2.47| ) into a diagonal form, let us consider the eigen- 
value problem for the non-differential part of ( |2.47| ): 

mm (2«W - v%> +1 - v^) = p(p + l)v%>, m = 1, • • • , JV-1, (2.48) 

where the eigenvalue is expressed as p(p + 1), and we have vjj^ = v$ = 0. The solution to 
eq. ( |2.48|) is given in appendix |B[ There we show that p characterizing the eigenvalue takes 
the integer values; p — 1,2, . . . ,N — 1. The explicit form of the eigenvectors v$ is also given 
there. Then, by expressing $ TO as 

N-l 

* TO (r)= E«SV W (r), (2.49) 
p =i 

the differential equations ( 2.47| ) are reduced to diagonal ones for <£?W; 



<pto\r)"- P r2 <P (p) (r) = 0. p=l,...,JV-l, (2.50) 
sinh r 

In order to solve (|2.50|) , we make a change of variables from r to y = e 2r , and consider 
{p( p \y) = (y — l) p ip^ instead of <£>^. Eq. (|2.50| ) is then transformed into a Gauss's hypergeo- 
metric differential equation: 

y(l - y)^ {p) + [1 - (1 - 2p)y] ^ (p) -pV (p) = 0. (2.51) 
The general solution to ( |2.50| ) is therefore given as a linear combination of two independent 



ones: 

M 



r ' =P P (y- l)- p (F(-p, -p, 1; ?/) + c p [F(-p, -p, 1; y) my + F*{-p, -p, 1; y)]) (2.52) 



where^] 



F(-P,-P,i;y) = E(^) w*. ( 2 - 53 ) 



fc=0 



F-(-p, -p, l; „) = -2 g Q ^ g (-j^ + I) (2.54) 

and /3 P and c p are constants. For our present purpose, c p must be chosen in such a way that 
(pW of ( P-52j ) is non-singular at y = 1 (r = 0). Such c p is given by 

-i 



/ p 1 < 

^2£ 7 j • (2.55) 



1 F* (a, /3, 7 ; y) = (8 /da + 0/8(3 + 2d/d~f) F(a, 0, 7; y) 
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The derivation of (|2.55| ) is given in appendix 0. There we obtain the solution using the 
analytic continuation of the hypergeometric functions. For example, ipW with p = 1,2 and 3 
are 



(i) 



j3\ (r cothr — 1) , 



(2) 



-/?2 coth r — r 



2 cosh r + 1 
3 sinh 2 r 



(2.56) 



(3) 



ft 



15 3 cosh r (2 cosh 2 r + 3) 

+ r— 



1 1 sinh 2 r 11 sinh 3 r / 

<y9 < - 1 ' ) and <y9^ 2 - ) are nothing but the solutions given in ]rjj as cp + and y_ for the AT = 3 case. 

As for the other scalar Y for the present (y a ) of ( |2.45| ), \l/ m of ( |2.33|) is given using ( |2.35| ) 
and ( ^46|) by 

^™ 

Since i;W oc mm, this is the p = 1 term of ( |2.49| ). 

3 Behavior of the solutions 



—mm (r coth r — 1) . 



(2.57) 



In this section, we shall discuss various aspects of the solutions: explicit analysis of the exact 
solutions for the cases N = 4 and 5 (sec. |3.1|) , the behavior of the solutions near the origin for 
general A" (sec. |3.2| ), the concept of the effective charges (sec. |3.3|) , and new solutions for the 
cases where some of y a are degenerate (sec. p^4] ). 

3.1 Exact solutions for 577(4) and 577(5) 

In this subsection, we shall study the exact solutions presented in the previous section in more 
detail for the cases N = 4 and 5. First, let us consider the 577(4) case. Since the eigenvectors 
v m' 2 ' 3 ^ ("with appropriate normalization) read (see appendix jB|) 



7( 2 ) 



1 





1 



7(3) 



(3.1) 



-1 



the X-coordinates of the four D3-branes are given by 







( 3 


1 1 \ 


x 2 


1 


1 


-1 -3 


x 3 


~ 2^ 


-1 


-1 3 


\xj 




\-3 


1 -l) 



(2) 
(3) 



(3.2) 
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in terms of f^ 1 ' 2 ' 3 ) of eq. (|2.56|) . The charges (e a ,g a ) and the D3-brane coordinates at infinity 
(x a ,y a ) determined from the asymptotic form of the present exact solution are 



(e 1 ,g 1 ) = (3(3 1 + (3 2 + p 3 , -3), (e 2 ,g 2 ) = (fa - p 2 - 3/3 3 , -1), 

(e 3 ,g 3 ) = (-/3 1 -(3 2 + 3/3 3 , 1), (e A ,g A ) = (-3(3, + (3 2 - (3 3 , 3), (3.3) 

and 




Figure 3: The trajectories of (X a , Y a ) of the N = 4 exact solution (solid lines). The parameters 
are chosen as (/3i, p 2 , j3 3 ) = (1, —2, —2), which give the electric and magnetic charges (e a , g a ) = 
(— 1, —3), (9, —1), (—5, 1), (—3, 3). The dotted lines describe the corresponding four-pronged 
string in the string picture. 

For (Pi, p 2 , p 3 ) = (~P, 0, 0), we have X a = PY a . In this special case, all the D3-branes and 
strings are located on a straight line which passes the origin. The charges (e a ,g a ) are parallel 
to each other and the configuration preserves 1/2 supersymmetries. 

In fig. 0, an example of the typical configurations is presented with solid curves. Notice 
that the four D3-brane surfaces bend and have a common tangent at the origin. We have 
chosen the parameters P p in order for the electric charges to take integer values. 
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4 : ,*i 4 ; 3 3 : ,*4 

• • • • • • 

• • • • • • 

• • • • • • 
* • • • • ^ *• * • • ^ * • ■ • • ^ 

• • • • • • 

• • * • • • 

3 • "'• 1 2 ""• 2 2 •' *• 2 

(a) (b) (c) 

Figure 4: Three different ways of connecting the four D3-branes. 



As explained in sec. \l.2[ we can find the configuration of a four-pronged string in the 
string picture from the asymptotic behavior of the solution. In the case of four D3-branes, 
there are three different ways to connect the D3-branes by a four-pronged string (see fig. 
Among them, the physical configurations are chosen by the condition that the tensions of all 
the strings should be positive. The four-pronged string in fig. || (dotted lines) corresponds to 
the case (a) of fig. f|. In general, the correspondence between the parameters j3 p and the type 
of the four-pronged string is as follows: 
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case (a) \/3 2 \ > 5a/^-|/%| or \f3 2 \ < — 1# 



11' 



? 3| 



case (b) 5^ ^\/3 3 \ < \(3 2 \ < 

case(c) ^\{3 3 \ < \fo\ < 5^|/9 3 |- (3.5) 

Note that the whole parameter space of /3 P is covered, namely, we can identify at least one 
physical four-pronged string configurations for every value of j3 p . An interesting fact is that 
in the regions 

two cases satisfy the positive tension condition. In the former region, the condition is satisfied 




by (a) and (b), and in the latter by (a) and (c). When \(3 2 \ = Sy^l/^! holds, the three cases 
degenerate. There is no internal string, and the four strings meet at a point, which is not 
generally at the origin. 

The analysis of the exact solution for SU(5) can be carried out in the same manner. In 
fig. H, we present the trajectories of the five D3-brane surfaces (X a , Y a ) and the corresponding 
five-pronged string configuration for a certain (3 p . Notice that the D3-brane surfaces 1, 2, 4 
and 5 have a common tangent at the origin (X, Y) = (0, 0), while the surface 3 sticks to the 
others. 
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Figure 5: The trajectories of the five D3-brane coordinates (X a , Y a ) of the exact solution (solid 
lines), with charges (e a , g a ) = (—1, —4), (—5, —2), (6, 0), (—3, 2), and (3, 4). The corresponding 
five-pronged string in the string picture is represented by dotted lines. 



3.2 D3-brane surfaces near the origin 

Since the two scalars X and Y must vanish simultaneously at r = due to the regularity 
of the solution, all the N D3-brane surfaces, which are described by the eigenvalues of the 
scalars, meet at the origin (X, Y) = (0, 0) of the transverse plane. In this subsection we shall 
study how the D3-brane surfaces meet at the origin for general N. 

For this purpose, let us consider the behaviors near the origin of the scalars X and Y in 
our exact solution of sec. |2.4| . Rather than expanding the expression ( |2.52| ) with respect to 
r, it is easier to return to the differential equation ( [2.50|) to know the behavior of ip^) near 
r = 0. In fact, eq. (|2.50| ) is approximated near r = as 



¥ i r ) ~2 V 3 ( r ) 







(r~0), 



(3.7) 



and the solution regular at r = should behave as ip^ ~ r p+1 . Therefore, the leading 
contribution to $ m ( p. 49 ) comes from the p = 1 term. This fact, together with the expression 
for v$ implies that $ m (r) ~ 2ax^nffir 2 {ax is an m-independent constant) and hence 



that the leading behavior of the eigenvalue of the scalar X is given by 



X m (r) ~a x (N -2m + l)r. 



(3i 
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The behavior of the other scalar Y should be the same as (|3.8|) except the constant ax since 
Y is also a solution to the differential equation ( |2.44|) . Therefore, we get 



(X m , Y m ) ~ (N - 2m + 1) (a x , a Y ) r (r ~ 0). (3.9) 



The concrete expression of the constant ay is obtained from eqs. ( |2.35 ) and the q m term of 



iV(iV 2 



T) Y,vl (3-10) 



From eq. ( |3.9| ) we can deduce the followings. The shape of the junction of iV D3-brane 
surfaces at the origin (X, Y) = (0, 0) of the transverse plane differs depending on whether N 
is even or odd. When N is even, all the N D3-brane surfaces have a common tangent at the 
origin: half of the D3-brane surfaces X m (1 < m < N/2) are "smoothly" connected to the 
other half (N/2 + 1 < m < N) of the surfaces at the origin (X, Y) = (0, 0). On the other hand, 
when N is odd, the above picture is true except for the D3-brane surface with m = (N + l)/2 
since, for this particular m, fl3.9| ) vanishes and the leading behavior of (X^+i)/2, Y(N+i)/2j is 
given by the next 0(r 2 ) term which is generically not parallel to (ax,C(y)- Therefore, among 
N D3-brane surfaces, N — 1 have a common tangent at the origin, while the remaining one 
meets with the others with an angle. We saw in sec. |3j] and in ref. fl(| that the above 



behaviors are actually realized in the cases with N = 3, 4 and 5 (see figs. ||| and [|). 

The above analysis about the behavior of the scalars near r = was concerning the exact 
solution of sec. |2.4| for X. To obtain the behavior of a general solution corresponding to 
(y a ) other than Q2.45|) , we have to analyze the differential equation ( 2.44j) near r = 0. Naive 



substitution of the leading expression of a 2 ^, a 2 ^ ~ mm, derived from ( 2.41|) into ( 2.44|) gives eq 



(|2.47 ) with sinh 2 r replaced by r 2 . Therefore, the same expansion as ( [2.490 leads to ( p.7|) and 



one would conclude that the above behaviors of the scalars near the origin are not restricted 
to the exact solutions. Though we have confirmed that this naive analysis is valid for N = 3, 
4 and 5, it must be checked that the other (p^ with q ^ p enters in a harmless way in the 
apparently higher order terms of eq. (|3.7|). 

Finally in this subsection, we shall give two comments. Our first comment is concerning 
the behaviors of the scalars near the origin when iV is odd: one D3-brane (X m , Y m ) with 
m = (N + l)/2 stick to the rest iV — 1 branes with an angle when N is odd. Noticing that 
the R-R charge g m vanishes only for this particular m (c.f. (|2.40| )), this phenomenon would 
be ascribed to the fact that we are describing the BPS states by the classical treatment of the 



SYM theory (i.e. electrically) as we mentioned in ref. [ 16[ . 
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Secondly, the fact that all the D3-brane surfaces have a common tangent at the origin 
when N is even means that all of them cannot be globally straight in the transverse plane. 
Otherwise, the charges (e a ,g a ) are all parallel. This is the case also when N is odd as is seen 
from the force balance condition. 



3.3 Effective charges 



From the discussions in the previous subsection, the solutions with non-parallel electric and 



magnetic charges are necessarily described by curved lines. As discussed in sec. |2]2| and 
appendix [A], the asymptotic behavior (r — > oo) of the curved lines reproduces the multi- 
pronged string configuration in the string picture. In this subsection, we shall point out that 
this feature can be extended to any finite r. 

The basic equations used in the proof of Appendix [A] arej] 



N 



E ( M a2/a - V a X a ) = 0, 



a=l 



{U a , V a ) 



{e a ,9a) 



U= \Q X + M Y \ =2tt 



N 



A' 



a=l 



0=1 



(3.11) 
(3.12) 
(3.13) 



Thus if we can define quantities at finite r which satisfy the equations similar to ( 3.111) , 
( 3.12|) and ( p. 13 ), then one can conclude the same results as in appendix |A|, that is, we can 
form a multi-pronged string in the string picture using the tangent vectors at r, and the 
configuration has the satisfactory properties concerning the tensions of the strings, the BPS 
bound, the formation of a planer tree junction diagram, etc. 

The natural generalization of x a and y a would be just given by the diagonal elements of the 
Higgs fields X a (r) and Y a (r), where we assume a gauge in which X(r) and Y(r) are diagonal. 



Viewing the asymptotic behaviors ( 2.11 ) and ( 2.12 ), the generalization of (u a ,v a ) would be 
-2r 2 (X^(r),F a '(r)), where the prime denotes the derivative in the radial direction. Then the 
generalization of (|3.11|) would be expected as 



W(r) = r 2 Tr (X'(r)Y(r) - X{r)Y'{r)) = 0. 



(3.14) 



It is enough to show ( |3.14 ) on the positive z-axis because of the spherical symmetry. Both 
X a and Y a satisfy the same differential equations ( |2.3U| ), which is in the form (rX a )" + 
J2h=i Aa(r){rXb) = with a symmetric matrix A a b (r). Using this fact and taking a derivative 
1 1 We take 9 = for simplicity without losing generality. 
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of ( |3.14|) , one finds W = 0. Since W = at r = oo from ( |3.11| ), the equation ( |3.14|) holds in 
the whole space. 

Now let us define the r-dependent charges (e a (r), g a (j)) by the electric and magnetic fields 
in the radial direction: 



diag (e a (r)) = 2r 2 £ r (r), diag {g a {r)) = 2r 2 B r {r) 



(3.15) 



where the diagonal form of the S r (r) and B r (r) is guaranteed by ( |2.25| ) coming from the 
spherical symmetry anzats. Then ( |3.12| ) trivially holds from the Bogomol'nyi equations ( |2.18| ) 
and ( p!9[) . 




Figure 6: The tangent lines to the D3-brane surfaces (solid curves) at r = oo (dashed lines) 
and at a finite r (dotted lines). 



As for the last equation (|3.13|) , we first define the energy in the spherical region \x\ < r in 
the D3-brane world volume by 

U r = [ SxU, (3.16) 

J\x\<r 



where U is the energy density appearing in fl2.2|) . Because of BPS saturation, the equality in 
(|2.4j ) holds for the energy U r and the charges Qx,Y{ r ) an d M^y^), the definitions of which 
follow (|2.3j ) for the spherical region \x\ < r. From ( |2.18|) and (|2.19|) , Qy(r) — Mx{r) = holds, 
and hence the two charges do not contribute to the BPS bound of U r . Plugging (|2.25|) and 
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the same equation for X into the definitions of the charges Qx{ r ) and My(r) and noting the 
fact that the second term of (]2.25|) does not contribute, they become 



Q x (r) = 47rr 2 Tr(X(r)X / (r)), M Y (r) = 47rr 2 Tr(F(r)r(r)). (3.17) 



Thus from the definition of the effective charges (|3.15 ) and the Bogomol'nyi equations (^T 



and ( |2.19|) , the same equation as ( p.!3|) holds for the quantities we have defined. 



Because of the regularity of the solutions at r = 0, it is a general property that the effective 
charges (e a {r), g a {r)) vary as functions of r and tend to zero as r — > 0. If the ratio e a (r) : g a {r) 
depends on r, the D3-brane trajectories bend. We show the SU(3) case in fig. |]as an example. 
The three straight lines starting at (X a (r), Y a {r)) in the direction — (e a (r), g a (r)) meets at a 
common point. Namely, the trajectories of the D3-branes on the transverse plane is developed 
in such a way that for every r the tangent lines to the trajectories form a string junction 
carrying the effective charges. 

3.4 Solutions with degenerate y a 

Up to now we have considered the cases with distinct y a . When some of the y a degenerate, we 
can construct solutions with magnetic charges different from ( |2.40| ) In this subsection, 
we shall discuss solutions with degenerate y a for the S77(4) case. 

Since we have y\ < y 2 < 2/3 < 2/4, the patterns of the degeneracy are classified as 

f(i) (s,s,s,s), 
(2/1,2/2,2/3,2/4) = < (ii) (-8,-8,-8, 3s) or (-3s,s,s,s), (3.18) 
[ (iii) Only one pair in (y a ) is degenerate, 

with s > 0. As stated in sec. |37|, the analysis of the solutions of the differential equations 



(|2.44j ) at r ~ does not depend on the values of y a . On the other hand, as we shall show 
below, the asymptotic behaviors of the solutions change drastically. 

Let us first examine the case (i). Since eq. (|2.37|) is ill-defined for degenerate y a , we 
introduce an infinitesimally small parameter 5 as follows: 

(2/1, 2/2, 2/3, 2/4) = (s + a ± 5, -s + a 2 5, s + a 3 5, s + a 4 6), (3.19) 



where J2t=i a a = 0. Substituting this into eq. (|2.37|) and taking the limit 5 — > 0, we obtain 
the asymptotic behaviors 

Qi(r)-^~ 2 re 2sr , Q^) ~ ^e is \ g 3 ( r ) „ -| re 2 -, (r - 00). (3.20) 
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It follows from these equations that the magnetic charges of this solution are given by (g a ) = 
(-2,-2,2,2). Then, substituting (|3~20D into (ggg) , the differential equation (fJ4|) in the 
asymptotic region is approximated as 

(2 -1 0\ 

.. = 0, (3.21) 

2 J 

where $ is a column vector <fr = $ 2 , $3) T , and the dotted part is multiplied by the powers 
of the exponentially decaying factor e~ 2sr . Since the eigenvalues of the 3x3 matrix in ( 3.21|) 
are and 2 (doubly degenerate), there are two independent divergent modes which behave 
as r 2 for large r. In order to obtain a physically sensible $ m , we have to adjust two of the 




parameters of the solution to eliminate the divergent modes. As we mentioned in sec. |3.2| , the 
number of free parameters in the SU(4) solution X regular at r = is three. Therefore, we 
are left with only one free parameter. Recalling that the other scalar Y is already a solution 
to the same differential equation as for X with regular behavior (|2.12j ) at r — > oo, we find that 
X must be proportional to Y, implying that the solution is an uninteresting one with parallel 
charge vectors. Analysis for the case (ii) leads to the same situation, namely, two divergent 
modes exist at infinity. Hence we conclude that for (i) and (ii) there is no non-trivial regular 
solution. 

Next, let us consider the case (iii). In this case, the magnetic charges of the solutions 
are (—2,-2,1,3), (—3,0,0,3) and (—3,-1,2,2) corresponding to three different degeneracy 
patterns of (y a ), y\ = y 2 , y 2 = Hz and y% = y^, respectively. We shall take the case y± = y 2 . 
Then, similarly to eq. ( |3.21| ), the differential equation ( |2.44| ) is approximated for large r as 

/2 -1 0\ 

* + ... = 0. (3.22) 
' \0 0/ 

In this case, the eigenvalues of the 3x3 matrix are (doubly degenerate) and 2, and there is 
only one harmful divergent mode. Therefore, contrary to the cases (i) and (ii), there can be 
non-trivial solutions since we still have two free parameters in the solution after eliminating 
one divergent mode. 

For these solutions regular at r — > oo we have 

(X 1 ,Y 1 )-(X 2 ,Y 2 )=0(\). (3.23) 



This is because eq. fl3.22j ) implies 



(2$i - $ 2 )" - ^ (2*1 - $ 2 ) + • • • = (3.24) 
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for 2$! — $ 2 = 2r (X\ — X 2 ), and the divergent r 2 term in the solution at r —>■ oo has 
been eliminated by fine-tuning to leave the 1/r term. From eq. (|3.23|) and the asymptotic 
expressions ( |2.11|) and ( |2.12|) , we find that (xi,yi) = (x 2 ,y 2 ) and (ei,gi) = (e 2 ,g2) hold for 
the present solution. In the IIB string picture, two strings 1 and 2 degenerate completely, and 
the situation is reduced to the three-pronged string. 

An example of the solutions for (y a ) = (—1/2, —1/2, 0, 1) obtained by a numerical method 
is presented in fig. [7[ We see that two D3-brane surfaces 1 and 2 approaches to each other as 
r — > 00 as mentioned above, though at the origin they behave in the way explained in sees. 



3.1 and K 



1 - 



. 75 





3 












Figure 7: A configuration with (y a ) = (—1/2,-1/2,0,1). The locations of two D3-brane 
surfaces 1 and 2 at r = 00 coincide. The behavior near r = is magnified in the small 
window. 



4 The existence condition of a three-pronged string and 
the force between two strings 

In the string picture, the existence of a three-pronged string connecting three parallel D3- 
branes depends on its two-form charges and the relative locations of the D3-branes [0. Let 



us consider three D3-branes denoted by Di, D 2 and D 3 , and a three-pronged string connecting 
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them (fig. H). We denote the two-form charges of the strings ending on the Di, D 2 and D 3 
by (ei,gi), (e 2 , g 2 ) and (e 3 ,g 3 ), respectively, where the charges are conserved at the junction: 
Yh=i e i — Yh=i 9i — 0. Since the force balance condition determines the relative directions 
of the three strings meeting at the junction, the shape of the triangle formed by the three 
D3-branes is constrained. One can show easily that the necessary and sufficient conditions for 
the existence of a three-pronged string connecting the D3-branes are given by the following 
conditions for each angle 9i at the vertices Di of the triangle: 

0i < Angle [{ej,gj), (e k , g k )} for different i,j, k. (4.1) 

Here Angle[i>i, v 2 ] denotes the angle between the two vectors V\ and v 2 having a common 
initial point. When the length of one of the strings vanishes, one of the equalities holds in 



( |4. 1| ) . In such cases the configuration is a stable BPS state with two strings located at the 
same point in the D3-brane world volume. When the D3-brane locations are changed further 
out of ( |4.1| ), the stability may be lost. Thus the process of changing D3-brane locations from 
( |4.1| ) to the outside may be observed as a decay process of a field theory BPS state of a 
three-pronged string to two dyon states |12L [l3|]. On the other hand, when two strings with 




Figure 8: The triangle formed by three parallel D3-branes (solid lines). The charge vectors 
are described by the lines with an arrow. The dotted lines denote a three-pronged string 
connecting the D3-branes. 



some two-form charges are located at a finite distance in the D3-brane world volume, and the 
configuration of the D3-branes satisfies ( |4.1| ), the strings may attract each other and decay to 
the stable BPS state of a three-pronged string. Thus the evaluation of the force between two 
dyons would be another way to see the conditions Q4.1|) . In this section, we will estimate the 



force between well-separated two dyons following essentially the approach of [[15], pT| , and will 
actually reproduce ([4.1|) . 
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Figure 9: Two strings with charges (e^gi) and (e 2 , g 2 ) (dotted lines) connect one common 
and two different parallel D3-branes (solid lines). The two strings are separated by R in the 
D3-brane world volume. 



Let us consider the following configurations of D3-branes and strings. Three parallel D3- 
branes D x , D 2 and D 3 are located at (xx,yi), (x 2 ,y 2 ) an d (x 3 ,y 3 ) in the transverse two- 
dimensional plane, respectively (fig. |9|). An (ei,gi) string connects D\ and D 3 , and an (e 2 ,g 2 ) 
string connects D 2 and D 3 . The two strings are separated by a distance R in the D3-brane 
world volume. The two strings appear as two dyons with the electric and magnetic charges 
(ei, <7i) and (e 2 ,g 2 ) under the different U(l) subgroups of SU(3), respectively, in the D3-brane 
world volume theory. The dyon corresponding to the (ei,gi) string is denoted by d\ and the 
other one by d 2 . We assume the distance R is large enough to allow treating the dyons as 
point-like particles. 

There act a Coulomb force between the two dyons. The Coulomb potential is obtained as 

71 

Vco U iomb(#) = ^(eie 2 + gm), (4.2) 

by simply adding the electric and magnetic contributions. 

Another force is generated by the long-range Higgs field. Let us first evaluate the mass of 
the dyon d\ by neglecting the effect of the other dyon d 2 . By calculating Qx,y and M X y from 
the asymptotic behaviors of the electric and magnetic fields (|2.11| )-( p.l4j) , and plugging them 



into fl2~4[), the mass of the dyon d\ is obtained as 



m l = 2irJe 2 1 +g 2 1 l 13 , (4.3) 



where /13 = J(x\— x 3 ) 2 + (yi — y 3 ) 2 is the distance between Di and D 3 . This mass agrees 



with the corresponding string mass by identifying 27ryef + g\ as the string tension. We shall 
evaluate the potential energy between the dyons by taking into account the independence of 
the distance /13 caused by the long-range effect of the other dyon d 2 . 
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From (|2 . 1 1| ) and ( |2.12|) , the long-range behaviors of the Higgs field generated by the dyon 
d 2 are given by 

X ~ diag(xi, x 2 , x 3 ) + ^ diag(0, u 2 , u 3 ), 

Y ~ diag(yi, y 2 , y 3 ) + dia g( > v 2, v 3 ), (4.4) 



where the u a and v a are related to the charges by ( [2.16| ). The i?-dependent terms of the third 



components cause the independence of the distance Z13. The factors are expressed explicitly 
as 



(u3,v 3 ) = \Jel + gl l 23 , (4.5) 
where l 23 is a unit vector parallel to the string connecting D 2 and D 3 : 

7 (x2~ x 3 ,y 2 -y 3 ) . , 

\j{x 2 -x 3 ) 2 + (2/2 — 2/s) 2 



Plugging these asymptotic behaviors into (|4.3| ) and taking the first order in 1/-R, we obtain 
the potential energy from the long-range Higgs field as 



7T 



V A Higgs(i?) = -^V e i + 9i V e l + 92 hs ■ I23. (4.7) 



Adding this to the Coulomb potential ([4.2|) , the total potential energy is obtained as 

V(R) = ~ {jel + glylel + gl l l3 ■ l 23 - ei e 2 - #i<? 2 ) • (4.8) 

The cancellation between the Coulomb and the Higgs contributions occur when the angle 
between l\ 3 and Z23 equals that between the two vectors (ei, gi) and (e 2 , §2)- In case the former 
angle is smaller than the latter, the force between the two dyons is attractive, and vice versa. 
Thus we obtain the inequality for 9 3 in ( [4.1[ ) from a field theoretical viewpoint. Repeating 
similar discussions for the other combinations of two strings, we obtain all the inequalities in 



5 Summary and discussions 

In this paper, we have constructed 1/4 BPS states in M = 4 SU(N) SYM theory by solv- 
ing explicitly the Bogomol'nyi equations under the assumption of spherical symmetry. The 
solutions correspond to the string theory BPS states of a multi-pronged string connecting 
iV different parallel D3-branes. Here each string ending on D3-branes carries the two-form 
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charges (ei, N— 1), (e 2 , iV — 3), • ■ ■ , (e^, — iV+1), where e a take real values satisfying J2 a e a — 0. 
The NS-NS charges e a (or electric charges in the SYM theory) may be quantized if the quan- 
tization of the collective modes of our solutions is performed. We have also shown that, by 
fine-tuning some of the parameters, we can construct solutions with magnetic charges different 
from above. 



As we have studied in sees. [2.2| , |3.3| and appendix |A], the trajectories of D3-branes have 



interesting behaviors. We have shown that, from the asymptotic behaviors of our solutions, 
we can generate the BPS saturated configuration of the corresponding multi-pronged string, 
that is, straight strings compose a tree diagram with junctions and the forces balance at each 
junction. On the other hand, the behaviors at finite r are quite different from the above IIB 
string picture. The trajectories bend non-trivially and all the D3-brane surfaces meet at the 
origin r = 0. Nonetheless, an interesting feature of the trajectories is that we can actually 
generate the string picture configuration at every r by defining appropriately the effective 
electric and magnetic charges. At the origin r = the directions of all the charge vectors 
with non- vanishing magnetic charges degenerate, and hence all the D3-brane trajectories with 
non- vanishing R-R two-form string charges meet there with a common tangent. This fea- 
ture apparently violates the SL(2, Z) duality symmetry, and may come from our technical 
preference that we treat the SYM theory merely classically, i.e. electrically |16|| . 

We have constructed exact solutions for general SU(N) in the case the vacuum expectation 
values of one of the Higgs fields (Y) are parallel to the magnetic charges. The number of the 
free parameters of our exact solutions are given by N, one from the overall factor of the Higgs 
field and N — 1 from the vacuum expectation values of another Higgs field (X) . In the general 
cases of our (non-exact) solutions, the vacuum expectation values of the first Higgs field can be 
varied continuously, and hence the total number of the free parameters is expected to be given 
by 2N — 2. The solutions should be regular at the whole space including r = and should 
converge to finite values at r — > oo. We have checked the number of the free parameters up 



to SU (6) by expanding explicitly the second-order differential equations ( |2.44|) at r ~ and 
checking numerically the safe convergence of the solutions at r — ► oo. A general proof of this 
fact has not been given. 

This number of the free parameters of our solutions is interesting because it does not agree 
with the expectation in the IIB string picture, as discussed in our previous work [IB|. Let us 



consider a tree diagram composed of straight string segments and three-string junctions in a 
two-dimensional plane. Here we assume the number of external strings ending on D3-branes 
is N and they have the same fixed R-R charges as the magnetic charges of our solutions. We 
have the freedom to take N — 1 NS-NS charges (the other one is determined by the charge 
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conservation). Then the relative directions of the strings in the diagram are determined by 
the force balance condition. We have the freedom to take the 2N — 3 lengths of each string 
segment in the diagram. We should not count the freedom of the parallel shift and the rotation 
of the whole diagram, because we fixed these degrees of freedom in obtaining our solutions 
(see section |2.3| ). Thus the total number of the degrees of freedom is given by 3iV — 4, which 
exceeds the above degrees of freedom of our solutions by iV — 2. We do not have any reliable 
explanation of this discrepancy, but may point out that the freedom of choosing the two-form 
(or electric and magnetic) charges does not seem to be fully incorporated as can be seen from 
the r ~ behaviors of our solutions discussed in the second paragraph of this section. Thus 
we hope the number of free parameters of such solutions may agree with the string picture 
if we find a better formulation in field theory respecting the SL(2, Z) duality symmetry to 
describe the 1/4 BPS states. 



As pointed out in section |3.1| , we found two different multi-pronged string configurations 
which can correspond to the same field theory BPS states in certain parameter regions of our 
solutions. This may question the one-to-one correspondence between the BPS states in string 
theory and those in field theory, and the resolution of this fact may be open. 

In section f|, we discussed the force between two strings connecting one common and 
two different parallel D3-branes by regarding the strings as dyons and calculating the long- 
range force between them in the SYM theory. In a certain parameter region of the D3-brane 
coordinates, the force is attractive, and a three-pronged string can be formed as a bound state 
of the two strings. This result is another support to the existence of a three-pronged string 
connecting D3-branes. 

The supermultiplet of a 1/4 BPS state should contain states with spins higher than one 
(see appendix |D|). The appearance of a field with such a high spin is unusual in a field theory 
without gravity. It is intriguing to note that a 1/4 BPS state is a non-local object and cannot 
appear as an elementary local field. To see this, let us consider a (0, 1) string connecting two 
parallel D3-branes separated by a distance I. This string state corresponds to a monopole in 
SU(2) SYM theory. One sees that the energy of a monopole is given by -E(o,i) = Tl/g, while 
its width is 5^(o,i) — (T/) _1 by viewing the explicit expression of the monopole configuration, 
where T and g are the string tension and the string coupling constant, respectively |[23|| . By 
performing the SL(2, Z) duality transformation on the equation i?(o,i)5r(o.i) = 1/g, we obtain 
E(p tq )5r(p tq ) = \p + qr\ 2 /lmT for a (p, q) string in general. Thus one cannot take 5r — > 
limit simultaneously for two strings with non-parallel NS-NS and R-R charges by taking an 
appropriate limit of r. This fact implies that a 1/4 BPS state cannot appear as an elementary 
local field in any SL(2, Z) equivalent description, because it has non-parallel electric and 
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magnetic charges. 

In appendix [E|, we have shown that our solutions are also solutions of the non-abelian BI 

| H a 



action proposed in [p2 |. Among the developments of the BI physics [|17], [2| 
remarkable fact is the coincidence between the brane worldvolume approach and the target 
space approach [p5|| . Since any solution of supergravity corresponding to string junctions has 
not been constructed yet, we hope our brane realization of string junctions will give a certain 
insight to a supergravity realization of string junctions. 

It is clear that there are some unsatisfactory points to be overcome in our formulation to 
describe the general 1/4 BPS states. We need more freedom to choose the magnetic charges. 

2?| , but how to incorporate loops in 



A BPS state of multi-pronged string can have loops 
our formulation is not clear. An SL(2,Z) invariant formulation is also needed. It is known 
that Nahm equations |28[ describe monopoles on D3-branes by 1-brane dynamics rather 
than the D3-brane dynamics as we did. It might shed light on the above problems if the 
formulation is extended to multi-pronged strings considering appropriate boundary conditions 
of 1-brane fields at junctions like those given in [30]. 

Since we use only two adjoint Higgs fields in our strategy of solving the Bogomol'nyi 
equations, our results may be also applicable to Af = 2 SQCD. Thus the states discussed in 
this paper might be related to the exotic states with non-parallel electric and magnetic charges 
observed recently in |3T| in M = 2 MQCD |32|] (See also |15|, |33| for some discussions on such 
exotic states). 

Note added: While this paper was in the final stage of preparation, an article |34| appeared 
which has an overlap with our discussions on exact solutions. 
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Appendix 



A String networks in the IIB picture 



In this appendix, we present a proof of the properties (A) and (B) of sec. ^]2| concerning the 
presence of the string networks for general N. The following proof is based on a substitution 
rule which allows us to reduce the number A" of the strings emerging from the D3-branes (we 
call these strings "external string" hereafter) by one. 

For any two among the A" external strings, say the strings 1 and 2, let us consider a new 
string 12 carrying the charge (ei 2 ,<7i2) = (ei + e 2 , g\ +g 2 ) ■ The transverse coordinate (£12,2/12) 
from which the new string 12 emerges is determined by the conditions that the BPS bound 
E'bps ( |2.4| ) of the energy and the angle 9 of eq. fl2.9| ) are common in both the original system of 



A" external strings (1, 2, 3, . . . , N) and the new system oi N — 1 external strings (12, 3, . . . , N). 
Then, this new string 12 satisfies the two properties: 

(I) The junction point of the strings 1 and 2 lies on the string 12. 

(II) Taking this junction point as the common end of the three strings 1, 2 and 12, the mass 
of the string 12 is equal to the sum of the masses of the strings 1 and 2. 

To show the properties (I) and (II), note first that the condition determining (#12,2/12) is 
expressed as 

c "(;::H'(;:)+ c *(;:)> ^ 

where C a (a = 1, 2, 12) is the 2x2 matrix given by 

C a =( 9a ~ 6a ). (A.2) 



e a g a 

Then, the property (I) claims that, for t-y and £2 determined by the first equality of 

*> V tr ("> ) = ( f 2 )+ M = (y )+ t u (y), (A.3) 



2/1 / V v 1 J \ 2/2 J \V 2 J V 2/12 / V ^12 

there exits ti2 which satisfies the second equality. Such t\2 is easily found by eliminating 
(x a , y a ) (a = 1, 2) from eq. flO]) using eq. (|A.3|) to get 



tuC 12 [ Uu )=t 1 C 1 (^)+t 2 C 2 ( U2 ), (A.4) 



V12J \ViJ V v 2 
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which is equivalent to 

*12 (e? 2 + £ 2 2 ) = *i (e? + 0?) + * 2 (4 + £2) • (A.5) 

This equation ( |A.5| ) also implies the property (II) since the tension and the length of the 
strings are given by T a = 2tt sgn(t a )y / e^ + g\ and i a = \t a \Je^ + g%, respectively: 

^12 = 7i4+T 2 £ 2 . (A.6) 

Now, the properties (I) and (II) implies that we can reduce the problem of proving (A) 
and (B) for A" external strings to that for A^ — 1 external strings. Since we know that (A) and 
(B) hold for A^ = 3, we conclude that they are valid for an arbitrary N. 



B Solving the eigenvalue problem ( |2.48|) 



In this appendix, we shall solve the eigenvalue problem ( 2.48|) , namely, we obtain the eigen- 
values p{jp + 1) and the corresponding eigenvector v$ . For this purpose let us introduce a 
polynomial V^'iz) of a variable z as 



JV-l 



(B.l) 



m=l 



Then, eq. Q2.48| ) is expressed as a second order differential equation for V^ p \z): 
ziz - 1) V (p) " + (z - 1) \N + 1 - (N - 3)z| V ip) ' 



+ 



N + l 



(iV - 1)2 -p(p+ 1) 



V® = 0, 



(B.2) 



where the prime denotes the differentiation with respect to z. In order to transform the 
differential equation flB.2j ) into a more familiar one, we rewrite it in terms of V^ p \z) defined 
by 

V {p) (z) = zil - zf-W^iz). (B.3) 
The differential equation for reads 



2(1 - z)V (p) " + [-N + 1 - (2p - N + 3)z]v {p) ' 

-(p- N + l)ip + l)V (p) = 0. 



(B.4) 



This hypergeometric differential equation has a polynomial solution when p takes integer values 
p = 1,2, . . . , N — 1, and we have 



V (p) (; 



z(l-z) 



N-p-l [N-p-l\ fk+p\ 

-l j2 v k M k I 



k=0 



(V) 



2*. 



(B.5) 
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The eigenvector v&> can be read off from (|B.1|) and ( |B.5|) . For example, v$ and vffl are given 
(up to the normalization constant) by 

v$ = mm, (B.6) 
vffl = mm (m — m) , (B-7) 

with m = N — m. 



C The coefficient c p 



In this appendix we derive the solution of (|2.51|) in an infinite Taylor series at y — 1 and 



obtain the coefficient c p in ( |2.55|) by an analytic continuation to the expression Q2.52|) . 

We should find the solutions (p^' of ( |2.51|) which damps faster than (1 — y) p at y ~ 1. 
Although the solution of the hypergeometric differential equation ( |2.51|) which is regular at 
y = 1 may be expressed formally by the hypergeometric series F(—p, —p, —2p\ 1 — y), this does 
not make sense because of a zero in the denominators. Thus the real solution may be obtained 
by taking an appropriate limit of ^ — > in the expression F(—p + e%, —p + e 2 , — 2p + e 3 ; 1 — y). 
The higher terms from order (1 — y) 2p+1 of this series have a 1/63 singularity at e3 = 0, 
and the series of these singular termsQ can be shown to satisfy by itself a hypergeometric 
differential equation which converges to ( |2.51| ) in the limit — > 0. Since this partial series 



linieg^o ^F{— p + €i, — p + e 2 , —2p + e 3 ; 1 —y) behaves like €162 at — > 0, we cancel this factor 
and obtain a solution 

{p {p) = lim lim —F{-p + e 1 ,-p + e 2 ,-2p + e 3 ;l-y), (C.l) 

ei,2— >0 €3^0 t\t 2 

where we wrote explicitly that the limit €3 — > should be taken first. 

A well-known formula for an analytic continuation of the hypergeometric function reads 
F(a, /3, 7; 1 - y) = r(7)r(7 - c*-/j) q + _ y) 

r(7-")r(7-p) 

+y^ a '" ni) l\\r^ l) F ^ - «, 7 - A 7 - « - /5 + 1; y). (C.2) 

We also use a well-known expansion formula for the gamma function: 

(-\) n (\ n 1 \ 
r(e -n) = ^J- (7 + E 7 + ^(1) + 0(e) J , (C.3) 



"Hence this series damps faster than (1 — y) p as desired. 
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where ^(1) is the Euler constant. Taking the limit in (|C.1|) by using (|C.2|) and ( |C.3| ), we 
obtain 



(p) 



(p! 



(2p)\ 



2 E 7 f (-p> ~v, i; y) + f*{-p, -p, i; 2/) + ln(y)F(-p, -p, l; y) , (C.4) 



W=l 



where we have used the definition F*(a, /3, 7; 2) = + + 2j^) F(a,f3, 7; 2). This is the 



solution ( gjg ) with c p = 1/(2 £f =1 



D Supersymmetric aspects of the solution 



It was argued in |L5| |1J| that the solution with non-parallel electric and magnetic charges has 
1/4 supersymmetry. In this appendix we shall examine this property in detail. 

The Lagrangian of the 3+1 dimensional M = 4 supersymmetric Yang-Mills theory is given 
by H 



-\g 2 ((fatcA^f + (f abc B\Bjf + 2 (f abc A»B; 



(D.l) 



where K and L are the indices of the 4 representation of the S77(4) .R-symmetry. The matrices 
a and /? span the antisymmetric part of 4 x 4, which transform in 6: 



0i 



0"! 



-CTl 



a 2 



0"3 



-0-3 



a 3 



V ^2 

-icr 2 



(D.2) 



Here <jj denote the Pauli matrices. The Lagrangian is invariant under the M = 4 supersym- 
metry transformation. The supersymmetry transformation of the gaugino is given by 



5\ 



K 



—I 



a^F^ - i 7 "D„ (a* A + ^TsSf 



+^S f -,ika k \A\ A 1 ] + ^ge,j t /} h [B', B>] + ga'/P [A', B'Yn,^ 



KL 



(D.3) 



where the parameter denotes a Majorana spinor. 



In investigating the supersymmetric aspects of our solutions, we identify (A3, .B3) with 
(X,Y) and put A\p and B\ 2 equal to zero due to the SO (6) (~ SU(A)) -R-symmetry of the 
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action. Plugging the Bogomol'nyi equations (|2.7|) , (|2.8| ), (|2.18| ) and ( |2.19|) into the above 
supersymmetry transformation (|D.3|) , we obtain 



5X1 



K 



ka 



(D.4) 



In our solutions the electric and magnetic fields behaves differently from each other, and thus 
we obtain two constraints^] on e^: 



KL 



KL 



<Ti 
°i 

01 

<Ti 



ia 2 



icr 2 J KL \-<Ji 



<n 



-io 2 
ia 2 



KL 





-a i 



e L = 0, 
e L = 0. 



Representing the Majorana spinor ti by a Weyl spinor t]l 

""'I VL ) ' 

and substituting (p. 61) into ( p.5|) , we find the conditions 



(D.5) 



(D.6) 



/ 1 i 
-i 1 

V 



1 -i 
i 1 / 



fVi\ 

\vJ 



(\ -% 

i 1 



The solution is obtained as 

Vl = a L rj, 



1 -i 
i 1 / 



(di, a 2 , a 3) a 4 ) = e ldR (0, 0, i, 1) 



fVi\ 

V2 
V3 

\vJ 



0. 



(D.7) 



(D.8) 



where rj is a constant Weyl spinor. Thus the BPS states discussed in this paper preserve one 
supersymmetry ( p.8|) out of the original four, and the phase factor in eq. (p.8| ) indicates the 
remaining U(l) i?-symmetry. 

In this M = 4 system, we have two independent central charges Qx and My after an 
appropriate £0(6) ^-transformation (a subgroup of this transformation is used to put 6 to 
zero in sec . The formula for the lower bound of the energy is given by 

U = Jd 3 x^Tr {{S^DiXf + (B lT D t Y) 2 + (D X) 2 + (D Y) 2 - [X, Y} 2 } ±Q X ±M Y 
> max{ \Q X + M Y \ , \Q X - M Y \ }, (D.9) 



^Our convention of the Dirac matrices are 
1 



7 



-1 



7 



75 



, ^ =2^,71- 
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which depends on the relative sign of the charges Qx and My. Though we solve only one of 
the two cases in this paper, this is not essential because the Bogomol'nyi equations of another 
case are just obtained by the substitution X — ► —X. Since our solutions have non-zero Qx 
and My, their masses saturate either | Qx + My| or | Qx — My|, and not both. Generally, 
states saturating k bounds are generated by 2{M — k) fermionic creation operators made of the 
supercharges. Hence our BPS configuration of SYM belongs to a supermultiplet with nx2 6 
(n: integer) components. This implies that this multiplet contains states with spins higher 
than or equal to 3/2. 



E The non-Abelian Born-Infeld action 



We shall show in this appendix that the BPS saturated solutions in this and our previous 



papers (16| are also solutions of the non-abelian Born-Infeld (BI) action proposed in |22] . We 
will follow the discussions given in |p3|| , where the usual monopole and dyon solutions in SU (2) 
SYM theory were shown to satisfy the equations of motion from the non-abelian BI action. 



The ordinary D3-brane effective action is given by the abelian BI action 

S = Jd A x ^/-det(^ + Fn 



(E.l) 



where we have ignored the overall numerical factor, which is irrelevant in the discussions here, 
and have appropriately normalized the gauge fields to be consistent with eq. ( |2.1| ). In our 
solutions, the D3-branes fluctuate in the transverse two dimensional plane with coordinates 
X and Y . Hence the induced metric on the D3-branes is given by 



g^u = Vnv + df,Xd v X + d,¥d v Y, 



(E.2) 



where t]^ denotes the flat Minkowski metric. The action (|E. 1| ) with the metric ( E.2 ) can be 
simply deduced [2M by the dimensional reduction from the 6-dimensional BI action with a flat 
metric: 



S = Jd 6 x ^-det(r] ab + F ab ) 



(E.3) 



Here the fluctuations of the four- dimensional induced metric ( |E.2[ ) are incorporated in the 
six-dimensional gauge fields in ( |E.3| ) by identifying X and Y with some components through 
T-duality: X = A 4 ,Y = A 5 . 

A non-abelian version of the BI action suffers from the ordering ambiguities of the non- 
commutative field components. This ambiguity is fixed in the proposal of |22| by taking the 
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symmetrized trace operation STr. The action is denned by 

S = J (fx STr^-det^ + Foj,), 

where the 6-dimensional field strengths are given by 
/ 



(E.4) 



ah 






£1 


£2 


£3 


D X 


D Y \ 


-£i 





B 3 


-B 2 


D X X 


D\Y 


-£2 


-B 3 





B x 


D 2 X 


D 2 Y 


-£3 


B 2 


-B x 





D 3 X 


D 3 Y 


-D X 


-D X X 


-D 2 X 


-D 3 X 





-i[X,Y] 


-D Y 


-DxY 


-D 2 Y 


-D 3 Y 


i[X,Y] 


/ 



(E.5) 



Since the part inside the symmetrized trace in eq. (|E.4f) can be expanded in a polynomial 
form in F a b, we can treat this part as if F a b were an abelian variable before symmetrization. 
The variation of the action gives 



5S=- STr 
2 



det^ + i^) ( V + F)- b l \5F 



(E.6) 



In evaluating ( |E.6| ) for our configurations, we are allowed to substitute eqs. (|2.7|), (^.8|), (P]T3) 
and Q2.19| ) into the curly bracket part {• • •} before the symmetrizing operation. This is because 
these four equations are linear relations among the components of F a b. Therefore, after a little 
algebra, the quantity inside the curly bracket of ( |E.6|) is simply replaced by — F a b. Hence we 
find a similar situation as ||23|1 , that is, there is no need to distinguish the symmetrized trace 
operation from the ordinary one, since we have only two non-commutative elements in the 
symmetrized trace of ( |E.6|) . Noticing that 5F a b = D[ a 5Ab] and performing an integration by 
parts, the equations of motion are reduced to 



D a F, 



ab 



0. 



(E.7) 



which are just the ones from the SYM action. The BPS conditions ( 2.7|) , ( |2.8|) , Q2.18 ) and 
(|2.19| ) satisfy eq. ( |E.7|) . Therefore we have shown that our BPS saturated configuration is also 
a solution of the equations of motion from the non-abelian BI action. 
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